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ABSTRACT 

In this article we consider the growth of seed black holes immersed in dark 
matter halos. We first investigate the adiabatic growth in various initial distribu- 
tion functions (isothermal, power law, and NFW) and find the resulting density, 
radial velocity, and anisotropy profiles. In addition we estimate the growth rate 
for a given black hole mass in the corresponding adiabatically modified dark mat- 
ter distribution function. Only in the isothermal case is there a convincing black 
hole mass-age relation. By calculating the line of sight velocity dispersion for the 
various cases as a function of the black hole mass, we find the predicted adia- 
batic M^h — c" relation; this never approaches the recently observed power law. 
We conclude by abandoning adiabaticity, suggesting that the black hole grows 
proportionally to the dark matter halo itself on a dynamic time scale. This allows 
us to relate the observed M^h — cr relation to the cosmological power spectrum 
on galactic scales by using dimensional scaling arguments. 

Subject headings: black hole physics — galaxies: evolution — galaxies: halos — 
galaxies: nuclei — dark matter 



1. INTRODUCTION 

Recent results by Gebhardt et al. (2000) and by Ferrarese & Merritt (2000) that establish 
a strong correlation between central black hole mass M^/j and the velocity dispersion a^-, 
measured at rg/S, indicate that the central black hole is intimately related to the dynamical 
structure of the galaxy. An earlier result by Magorrian et al. (1998) relates Mf,/i linearly 
to the mass of the bulge, which suggests a similar conclusion, while Merrifield, Forbes, & 
Terlevich (2000) establish a link between M^h and the age of the stellar system: a massive 
central black hole seems to grow over a period of a Gyr or so. 

These results have already stimulated the emission of various theories that create a 
feed-back mechanism between bulge star formation and accretion of gas onto the black hole 
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(Burkert & Silk 2001). These theories generally tend to estabhsh the desired correlations 
but at the cost of some rather complicated physics simply, perhaps oversimply, described. 

The object of the present paper is to first revisit the correlations estabhshed by the 
adiabatic growth of a black hole in a galaxy, since this process is relatively free of physical 
assumptions once the initial distribution function (DF) is chosen. This part of the work is 
very much in the spirit of van dcr Marcl (1999), who examined the effect of the black hole 
on the fundamental plane relations, but did not consider explicitly the — a relation. 
We choose distribution functions moreover that are appropriate for collisionless dark matter 
halos, the philosophy being that these are the dynamically dominant components of massive 
galaxies and should therefore dictate the observed velocities. These are chosen to be an 
isothermal or Gaussian DF (to test our code against previous work and because this may be 
the maximum entropy state according to Nakamura (2000)), an isotropic steady-state power 
law DF found by Henriksen & Widrow (1995) for comparison purposes, and the DF that 
corresponds to the NFW density profile, as approximated by Widrow (2000), as the best 
measured approximation to a dark-matter halo. 

Such an approach does not, of course, explain the origin of the black hole, but rather 
yields only the perturbed DF that is created by its adiabatically estabhshed presence. An 

earlier attempt to grow the black hole during the formation of the galaxy (Ostriker 2000) 
relied on dissipative dark matter which is fraught unfortunately with badly known parame- 
ters. In the present approach the black hole growth is limited to the particle fiux across the 
event horizon of a seed black hole that is peculiar to the initial or adiabatically modified DF. 

The adiabatic approach fails to yield either the correct form or extent of the observed 
correlations between the black hole mass and the modified galaxy. We therefore suggest a 
possible explanation of the correlations based on a non-adiabatic process of black hole growth 
on the formation timescale of the dark matter halo. The argument is essentially dimensional 
at this stage and must be examined numerically in greater detail. 

In section 2 we review the adiabatic growth approximation and verify our code with 
the isothermal DF. In section 3 we give the results in the power law and NFW distribution 
functions. In section 4 we discuss the resulting M^/j versus in the various distribution 
functions and show the extent of the black hole influence in the galaxy. In section 5 we 
give our dimensional derivation of the observed relations based on a self-similar (but not 
necessarily spherically symmetric) growth of black hole and dark halo. Finally we give our 
conclusions. 
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2. ADIABATIC GROWTH 
2.1. History 

We consider the slow growth of an initially small seed black hole (BH) located at the 
centre of a coUisionless, spherical system of stars and dark matter. Using an algorithm that 
makes use of the fact that, under these circumstances, the radial and transverse actions are 
conserved, the final state of the system containing a supermassive BH can be calculated. 

This technique was first suggested by Peebles (1972), who used it to show that an 
adiabatic cusp with q ~ would form in an isothermal sphere. This work was confirmed 

numerically by Young (1980). His algorithm is used here essentially unchanged. Lee & 
Goodman (1989) used Young's algorithm to explore the more general case of a stellar system 
with a net rotation, although they treat the potential as spherically symmetric rather than 
axisymmetric. They found that, while the rotation to dispersion ratio v/a was larger after 
adiabatic growth, the gain was not enough to match the observed values. A variety of 
systems were examined in detail by Quinlan et al. (1995), and van der Marel (1999) used 
their code to explain a number of observational properties of galaxies. Finally, the adiabatic 
growth model was used by Gondolo & Silk (1999) to examine how dark matter annhilations 
in the cusp that would form in our own galaxy could be used as a probe of the dynamics 
or of the nature of dark matter. They suggest that searching for a neutrino signal from the 
central density spike (the neutrino flux increases with an increase in cusp slope) could set 
upper bounds to the dark matter cusp in our galaxy, or, alternatively, clarify the nature of 
dark matter. 

Although these studies serve to demonstrate the utility and topicality of the technique, 
they have not studied explicitly the relation between dark matter distribution functions and 
the central black hole mass that is our concern below. 
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2.2. Units 

For convenience, we use here dimensionless units, characterized by the variables 

R = r/ro, 

Q = P/Po, 

e = iE-%)/iA7rGporl), 

= (<l>-<l>o)/(47rGpor^), (1) 

m = M/(47rporg), 

The fiducial quantities Tq and po are usually taken to be the core radius of the galaxy, similar 
to the usual definition of the King radius, and the density at this radius (where the surface 
density is about one half of its central value in a fitted isothermal sphere). 



2.3. Theory 

The initial system, containing the negligible mass BH, can be described by a distribution 
function F{e,j), and the density of the system can be calculated from 

g{R)^A7r de ^^^^^i), (2) 

Jo ^ 

where the radial velocity is expressed as Vr = \2{e — '4>) — j^/i?^]^/^ and j^ax = [2R'^{e — ^jJ)Y^'^ 
is the maximum angular momentum. The potential of the system can then be calculated 
from Poisson's equation, 

^d^rdsj=^- 

The integral form of Poisson's equation is often more convenient, however. We take the 
reference potential to be zero at the centre of the system and we suppose the halo mass to tend 
to zero there, which is true for a density profile less steep than R~^. Hence ip — dip/dR — 
for i? = 0, and we have 

^{R)= £g{s)sds-^, (4) 

where the total mass of the system is 

m{R) = / g{s)s'^ds + rribh, (5) 
Jo 
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and rribh is the mass of a central object, if one exists. 

Now, if the growth of the central seed black hole is much slower than the dynamical 
time of the system, the radial and transverse actions of the system are conserved. These 
so-called adiabatic invariants are given by (Young 1980) 



iR = 2 VRdR 

where R± are the turning points of the orbit (zeros of the integrand) , and 

^T = CVrdiRO) 
= 27rj, 

so that the angular momentum is an action variable. 



(6) 



2tt 

(7) 



The adiabatic growth framework is then as follows. A particle of the initial system with 
energy e and angular momentum j will have a different energy e* (but identical angular 
momentum j) after the BH has grown. However, its radial action ir will be the same before 
and after growth; thus, it is possible to invert i\{e* ,j) to get £ as a function of (£*, j). The 
final DF is then found from F*{e*,j) = F{s{s*,j),j). 

This is implemented in an algorithm similar to that used first by Young (1980). See 
Appendix A for details. 



2.4. The Isothermal Sphere 

We apply the adiabatic growth framework first to the isothermal sphere, characterized 
by the Gaussian DF 

F(£) = (27r)-3/2g-e_ 

This is a well-studied distribution, and was Young's (1980) initial system. We present 
the adiabatic results for the isothermal sphere here mainly as a check that the code devel- 
oped gave the same results that Young (1980) found. There are in fact reasons to beheve 
(Nakamura 2000; Henriksen & Le DeUiou 2001) that this is ultimately the DF of interest, 
especially near the centre of the system. Of course the system may not be fully relaxed when 
the black hole is growing. 

Results are shown in Figure 1 for black hole masses from zero through 0.001, 0.01, 0.1 
and 1.0 as labelled by increasing effect. These dimensionless masses correspond to physical 
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masses in the range of about 10^ to 10^° Mq if we use the values of po and vq given in section 
3.3. 

We find the same characteristic density cusp of as obtained by Young (1980), and 

the radial velocity profiles are also similar. Although Young (1980) did not show results for 
the anisotropy parameter, defined as P — 1— < > / <2T/|>, our results do agree with 
the prediction of Quinlan et al. (1995) for distributions of this type. Wc observe that there 
is at most about a 10% anisotropy in favour of tangential motion, due to the increasing 
binding energy of a particle and the resultant decrease in eccentricity at constant angular 
momentum. Moreover the perturbation extends as far as the core radius only when the mass 
is comparable to that of the core, as is to be expected. These results are the main test of 
our program. 

3. CUSPY DARK MATTER HALOS 

We explore here the adiabatic growth of a central black hole in a collisionless dark 
matter halo that possesses a central cusp. For initial systems, we choose two very different 
starting points. The first is an isotropic self-similar system, meant to represent the final 
state of a halo undergoing self-similar relaxation (Henriksen & Widrow 1995). The other 
system is the NFW system of Navarro, Prenk & White (1996), which fits a wide range of 
dark matter halo sizes. 



3.1. Self-similcir Distribution 

The self-similar system is described by the DF 

F(£) = Fo|£|"i^, (9) 

where 5 is a free parameter (2/3 < 5 < 1, 5 > 1) that essentially controls the logarithmic 
slope of the initial density and potential, given by 

Q = R-^l^ and (10) 

^P = (11) 

The constant Fq depends only on the value of 5, and can be solved for in terms of the core 
radius Tq and the density po by using equations (2) and (4). 



- 7- 



Using this system as a starting point for the adiabatic algorithm discussed above leads 
to rather different results than are found for the isothermal sphere. The results for two 
different choices of 6 are given; Figure 2 is for the parameter 5 — 2, while Figure 3 shows 
results with 6 = 3/4. 

As can be in Figure 2, the initial density profile of the self-similar system with 6 = 2 
is a power law with a logarithmic slope of —1. Adding an adiabatically grown BH to the 
system induces a steeper cusp region at small radii, and the outer radius at which this region 
begins depends of course on the final BH mass. The new cusp slope is —7/3, greater than 
the density cusp induced in the isothermal sphere, but the difference between the bcforc- 
and- after slopes are not as great here as in the isothermal case (a gain of ~ 1.3 here versus 
1.5 for the isothermal sphere). The major difference with the isothermal sphere however is 
best seen in the anisotropy parameter. The system becomes quite tangentially anisotropic 
at small radii. 

That the cusp slope will not increase as greatly for an initial system which already has 
a steep density cusp is shown best in Figure 3, where the grown BH seems to have left no 
visible mark on the density slope (to within graphical resolution: the initial density slope is 
—8/3 while the measured final slope is —11/4). This result is for 6 = 3/4. Notice, however, 
that the BH growth has disturbed the radial velocity, which takes on the same 1/R shape 
(although only for the largest masses) seen previously in both the 6 = 2 self-similar system 
and the isothermal sphere. That these systems all have similar velocities at small radii is 
simply indicative of their similar Keplerian potentials there. 



3.2. The NFW Profile 

The NFW system, given by the universal density profile 

is more complicated to describe than the other two systems discussed above, since it does not 
have an analytic distribution function. We use here an analytic fit of a numerical calculation 
of the DF (Widrow 2000), which in our units is 



F{e) = F,{1 - ef'e-'/' (^-^1^^ e^. 



(13) 
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where P — ^^PiS^ and the parameters are: 

Fi = 9.1967 X 10-2 
q = -2.7419 
Pi = 0.3625 

P2 = -0.5669 ^ ^ 

P3 = -0.0802 
P4 = -0.4945. 

The results for this system are given in Figure 4. As expected, the BH growth induces 
cusps in both the density and velocity. The logarithmic slope for the density is —7/3, which 
is the same as for the self-similar system with S — 2; they also both began with the same 
slope oc l/R) in the inner region where the BH disturbance is greatest. The velocity cusp 
is the usual l/R. We notice that the velocity tangential anisotropy is also comparable to the 
system with S = 2 at small radii for all black hole masses, and that the velocity perturbation 
can extend nearly to the core radius for rribh > 0.1. 



3.3. Growth Timescales 

It is of some interest to explore the rate at which a black hole would grow in the 
above dark matter distributions in view of the black hole mass galactic age relation reported 
by Merrifield et al. (2000). Such growth represents an alternative to advection dominated 
accretion flow (ADAF; see Narayan, Igumenshchev & Abramowicz (2000) and references 
therein) as a means to grow a black hole invisibly. We consider as an illustration the time 
to grow the black hole by a factor of ten in the various DFs. 

In order to calculate the timescale for the BH to grow by a factor of ten, we use the 
following expression for the growth of the BH as it accretes matter: 

= ^TTrlp{r^) <Vr> \r,, (15) 

where Vx is the radius of the last stable orbit, = SRg, and Rg — 2GMbh/c^ is the 
Schwarschild radius of the BH. In our units, and with < Vr > given as an integral over 
the distribution function, this expression becomes 

^ = STT^AnGy/'pfrl / de / jdjF{e,j), (16) 

JrP, Jo 

where = "^{fx) is the potential evaluated at the radius of the last stable orbit. Note that 
this equation requires the central density po and the core radius Tq. These two parameters 
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control the size and shape of the system; for simphcity, we choose the values 

To = 500 pc, 



Po = 30Mepc-3; ^^^^ 

these values are averages of the systems van der Marel (1999) used in fitting the adiabatic 
growth model to a variety of different galaxies. 

We carry out the calculation (16) numerically, taking the DF dependence on black hole 
mass to be given by assuming adiabatic growth. This was done for 30 different black hole 
masses, resulting in an equation of the form 

= i^Mj^h, (18) 

with b close to 2 for the isothermal system, but rather different in the other cases. Such an 
equation is readily integrated to give the growth time scale. The mass is not subtracted from 
the DF itself as both our trials and those of Quinlan et al. (1995) suggest that the effect on 
the results is small, at least for rubh < 0.5. 

We first report the calculation for the isothermal sphere (containing an adiabatically 
grown BH). Numerical calculation of equation (16) gives b — 2 and k, — 10~^^ Mg^yrs"^; 
integration of (18) with a seed BH mass of IO^Mq then results in a timescale of about nine 
billion years. This is close to what is interesting, since this is about the time available since 
the quasar epoch. 

If we instead use the self-similar DF, modified by the adiabatic growth of a central BH, 
the timescale becomes absurdly short. For 6 — 3/4:, we calculate 6 ~ —1/2 and k, — 10^^ 
Mq'^-Yts''^, leading to a very rapid growth (using again a seed mass of 10^ Mq): it would 
take only 10"^*^ years to grow the BH by a factor of ten! The DF for 5 < 1 increases with 
increasingly negative energy and the density profile is very steep. Moreover, the orbits tend 
to remain rather isotropic as the BH grows, so that essentially all of the mass falls in on a 
central dynamical timescale. Setting 6 > 1 does slow this growth down substantially; for 
5 = 2, for example, we get k = 250 Mg^^yrs"^ with b = 1/3, giving a timescale of about 
5,000 years. Increasing 5 increases this time to at most a few hundred thousand years. 

The NFW system (modified by the BH growth) is similar in some respects to the self- 
similar system with 6 = 2; both share the same density and radial velocity profiles at small 
radii. Indeed, the timescale calculation for the NFW system yields results similar to the 6 = 2 
case, with the only difference being a smaller value for constant, at k, = 115 Mq^yts"^. This 
smaller constant gives a timescale of just over 10,000 years for the BH to grow ten times 
larger than its initial mass of 10^ Mq. This may be accounted for by the more pronounced 
anisotropy of the NFW DF at larger radii. 
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It should be noted that changing the initial constants (ro, po, and the initial BH mass 
Mq) can vary these numbers by a few orders of magnitude in both directions. From the 
work of van der Marel (1999), which fits an isothermal sphere to a variety of galaxies, it 
seems that the core radius can be as small as a few tens of parsecs for "power-law" galaxies 
and as large as a few kiloparsecs for "core" galaxies. Assuming the same scaling laws van 
der Marel (1999) uses, the density po will depend on the value we take for the core radius. 
According to these scaing laws, a radius of Tq = 20 pc, for example, corresponds to a denisty 
of po — 8360 Mq pc~^, while at Tq = 2000 pc the density will he po — 3 Mq pc~^. Using 
these two sets of values, and assuming that the initial seed BH mass can vary between 10^ 
Mq and 10^ Mq, we see that the isothermal sphere can grow by a factor of ten anywhere 
between a million years and 10^^ years. For the self-similar system with S < 1, changing the 
constants has essentially no effect, since the timescale is too short; for 5^1, however, the 
timescale can approach a few billion years. The timescale for the NFW system will not go 
beyond a few hundred million years. 

It seems, then, that only the isothermal or Gaussian DF can yield a black hole mass 
age relation of the type detected by Merrifield et al. (2000) by adiabatic growth from a 
collisionless DF. The other systems studied here grow more quickly than the isothermal 
sphere, mainly because of their stronger central density cusp. 



4. THE ADIABATIC M^n - o RELATIONSHIP 

Recent observations (Ferrarese & Merritt 2000; Gebhardt et al. 2000) have found a 
strong correlation between the mass of the central BH and the line-of-sight velocity dispersion 
in the bulge of its host galaxy. Although various theories have been suggested to explain 
this relationship (e.g., Haehnelt & Kauffmann (2000); Adams et al. (2001)), none have yet 
to be proven conclusively. This relationship has been shown to follow 

M,h oc a", (19) 

where a is somewhere between 3.5 and 5. 

The adiabatic growth model has been used by van der Marel (1999) to explain various 
observational properties of black holes, such as the central density cusp and its correlation 
with the luminosity of the galactic bulge. His analysis included properties both intrinsic to 
the adiabatic growth as well as scaling relations based on fundamental-plane-like observa- 
tions. We repeat his analysis here to explore the M^h — cr relation. 

First we calculate any intrinsic relation between the BH mass and the velocity dispersion 
of the bulge that may arise naturally from the adiabatic growth of the central BH. Keep in 
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mind, however, that the calculations that follow cannot be compared directly with observa- 
tions, since they are noise-free and have an infinite resolution. Regardless, they should give 
a sense of whether or not the BH growth can give a relation like (19). 

The line-of-sight velocity dispersion is found by projecting the radial and transverse 
velocity moments on the plane of the sky. The velocity moments are calculated from 



< V^V^ >-- 



An 



Q J^iR) Jo ^ ^R 



projecting them on the sky gives 

2 /-^ dRRp 



1 _ s 

R^ 



<vl> < > 



2i?2 



^R,) Jr^ 

where Rp is the projected radius, and E is the projected density, given by 

dRRp 



JR 



(20) 



(21) 



(22) 



It is a simple matter to predict the dispersion near the centre of the system. As stated 
above, the velocity moments simply reflect the Keplerian potential near the BH; thus they 
take the form 



< V >cx 



R ' 



(23) 



and this form is identical regardless of the intial system. Writing the dispersion as cr = [< 
y2 >]i/2^ ^jijg relation is simply 



TJibh = o- . 



(24) 



Of course, this relation is applicable only in the innermost regions. Observations are 
usually done much farther from the centre - typically near the effective, or half-light, radius 
of the bulge. This radius corresponds to the core radius Tq for the unperturbed isothermal 
sphere. To calculate any relationship between the BH mass and the velocity dispersion away 
from the centre of the system, we must use the DF calculated from the adiabatic growth 
framework. 

Figure 5 shows the results for the three systems studied, calculated at three different 
radii: R — 10~^, near the centre; i? = 1, at the core radius; and R — 100, well outside of 
the central region. The simple predictions made above (24) are confirmed to exist in this 
model asymptotically as the mass of the hole becomes large (the solid fines in Figure 5). As 
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the calculations are done farther out, however, the relation exists only at larger and larger 
masses (the dashed and dotted hues in Figure 5). 

Similar calculations for the self-similar and NFW distributions yield almost identical 
results. The NFW DF is notable for attaining the asymptotic relation at smaller black hole 
masses than for the isothermal case (the self-similar DF is intermediate), but in all cases 
no linear relation steeper than that of equation (24) is found to exist. There is, as can be 
seen in the figures, a steep shoulder during the approach to the asymptotic limit but this is 
non- linear and of insignificant extent. 

So it seems that the required relation does not arise naturally in the adiabatic growth 
framework. This is not suprising, of course; after all, the BH does not create a disturbance 
much further than the radius R = rnhfi. Thus only the largest black holes can reach out to 
the typical radius at which observations are taken, and then they establish a much flatter 
relation. 

In the next section therefore we consider an alternative to adiabatic growth wherein the 
central black hole and the dark matter halo form together on the dynamical time scale. 

5. DYNAMICALLY GROWN BLACK HOLES 

The argument in this section is somewhat more speculative than in the preceding sec- 
tions, and it must ultimately be checked by extensive numerical calculations. However the 
argument is compelling on dimensional grounds and is consistent with well known sohitions 
and simulations in spherical symmetry (Henriksen & Widrow 1999; Fillmore & Goldreich 
1984; Bertschinger 1985). 

We assume that the galaxy forms by the extended collapse of a "halo" composed of 
collisionless matter, and that simultaneously the central black hole is growing proportionally 
to the halo as matter continues to fall in. We do not assume spherical symmetry. 

The preceding assumption is equivalent to the assumption of multi-dimensional self- 
similarity as defined in (Carter & Henriksen 1991) and in Henriksen (1997). The technique 
was used in spherical symmetry in Henriksen & Widrow (1995). The essential idea is that 
under this assumption one can write the mass inside any surface (the surface in space would 
be defined by holding M constant at a fixed time) in the halo as 

M = M(X)e(='^-'")^. (25) 
Here the vector X = re~^^'^^ is a scaled position vector and e*^"^) = at. The quantities 
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5 and a are scales in space and time that allow dimensional information to be included in 
the expressions for the various quantities. Thus the mass scale n is determined in terms of 
S and a by the condition that G is a constant of the problem. This yields /i — 3S — 2a as 
used above. 

The velocity of any particle in the halo may be written consistently as 



The quantities X and Y are independent of T during the self-similar collapse and so 
define a steady-state phase space. Consequently on a fixed spatial surface and averaged over 
the line of sight rms velocity we can eliminate T between equation (25) and the averaged 
equation (26) to obtain 



where a is the velocity dispersion along the line of sight. It does not matter which mass 
surface is chosen in the system if it is truly self-similar of course. Thus the preceding relation 
applies to the "bulge" mass at the bulge scale and to the black hole mass on the black hole 
scale. However the black hole mass will be simply proportional to the bulge mass during self- 
similar growth so that both masses will obey relation (27). This of course also requires that 
the black hole mass and the bulge mass are also proportional but, given that we are talking 
about total masses and recalling the vagaries of star formation, this would not necessarily 
imply a tight bulge luminosity black hole mass correlation. 

We may proceed to require that the constant of proportionality in equation (27) is equal 
to the observed (Gebhardt et al. 2000; Ferrarese & Merritt 2000) constant, say a. Then we 
find that 



so that, for example, if a = 4, then 8/a = 2; and if a = 4.5 then 5/a — 5/3. Should a < 4 
(but > 3 which appears here as a kind of lower permissible limit) say 15/4, then 5/a — 7/3. 
The principal numerical fact to note is that 5 /a is greater than or less than 2 depending on 
whether a is less than or greater than 4. 

The reason for the numerical discussion of the preceding paragraph is that the value 
5/a — 2\s highly significant in spherical models of dark matter halo growth (see e.g. Hen- 
riksen & Widrow (1999) for a summary). In these models this ratio is given in terms of the 
power law index — e of the initial cosmological density perturbation by 



V = 



y^g(<5-«)T_ 



(26) 




(27) 
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Consequently we can infer from equations (28,29) that the initial cosmological overdensity 
had the power — e where 

e^2(l-i). (30) 

If finally we relate e to the power spectrum index n of the primordial density through 
the rms profile of such perturbations (other choices are possible, e.g. Hoffman & Shaham 
(1985), but similar results are found), then 

n = 2e - 3. (31) 

Consequently we arrive at a direct link between a and the primordial power spectrum index 
on the scale of galaxy halos as 

n=l . (32 

a 

Thus under our interpretation of the black hole mass- velocity-dispersion correlation we 

are led to conclude that the primordial power spectrum on the scale of galaxies has the 
index n of equation (32). For a — A this yields n — —2, while a — 15/4 and a = 9/2 yield 
n = —11/5 and n = —5/3 respectively. These values for n on the scale of galaxies are all 
in good agreement with observation, which favours a value near n = —2. We conclude that 
our interpretation of the mass-velocity-dispersion relation as originating in the primordial 
density profile is consistent with cosmological evidence. 



6. CONCLUSIONS 

In this paper we have explored the implications of growing a black hole in various dark 
matter halo distribution functions. Our principal approach was to assume that the black 
hole grows adiabatically on a time scale long compared with the dynamical time of the halo. 
The method gives definite predictions for the modified density, radial velocity and anisotropy 
profiles in the isothermal, self-similar "power law" and NFW dark matter distributions. The 
isothermal calculations reproduced and extended slightly previous work, but the calculations 
for the self-similar and NFW dark matter halos are new. Depending on the mass of the black 
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hole the disturbances can be noticeable out to nearly the core radius of the galaxy. Moreover, 
estimates of black hole growth time scales in the adiabatically modified distribution functions 
arc given. Only in the isothermal DF is there found a reasonable black hole mass-galactic 
age relation. In no case however does the adiabatic argument give an Mj,h — a relation that 
is close to that observed. 

Thus in the concluding section wc explored by dimensional argument based on the 
concept of multidimensional self-similarity (Carter & Henriksen 1991; Hcnrikscn 1997) the 
possibility that the central black hole grew on a dynamical time scale with the dark matter 
halo. A simple argument predicts a = 3/(1 — e/2) where a is the observed power in the 

— o relation and — e is the power in r of the initial cosmolgical density perturbation that 
produced the galactic halo. Under certain assumptions this power can in turn be related to 
the power n of the primordial cosmological power spectrum on the scale of galactic halos, so 
that reasonably, a — 12/(1 — n). This gives n = — 2 for a = 4, which is close to that observed 
in both cases. We conclude that this suggestion is promising and more work should be done 
to confirm or infirm the idea that the black hole can grow self-similarly with the dark matter 
core. 

This work was supported by the Natural Science and Engineering Research Council of 
Canada. We thank the anonymous referee for helpful comments. 



This appendix describes the framework that was used to calculate the adiabatic growth 
models presented above. 

We implement this framework using an algorithm similar to Young's (1980), which solves 
directly for the final state of the system containing the supermassive black hole. Given a DF 
that describes the intial system, we 

1. Compute the self-consistent potential and density q for the initial system (eqs. (2) 



A. Algorithm 



and (4)). 



2 



Compute the radial action iji{e,j) (eq. (6)) for the intial potential ip. 



3 



Approximate the final potential for a given black hole mass by 



rribh 



(Al) 



R 
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4. Calculate the radial action i*j^{e,j) for this new potential ip*. 

5. Equate the intial radial action iji{e,j) and the final action i*ji{e,j) to find the energy 
e that has become e*. Then the new DF will be 

F*{e*J)^F{e{e*J)J). (A2) 

6. Compute a self-consistent density g* and new potential iJj* using F*. Continue back 
to step 4 and iterate until the density has converged. 

The convergence criterion for the algorithm was taken to be when g* changed by less 
than 10~^ at all radii. 

The computer program developed to implement this algorithm, written in C, uses a 
radial grid Ri, with points spaced logarithmically between an inner radius of i? = 10~^ 
and an outer of i? = 10^. The various properties of the system are described on this grid; 
gi = g{Ri), ipi = ip{Ri)-i and so on. The DF is described on a grid of energy and angular 
momentum points, where the energy points have for convenience the values of the potential 
on the radial grid, £j = ■^j, and the angular momentum is evaluated &s x — j/ jc, where jc is 
the circular angular momentum and x goes from zero to one. 
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Fig. 1. — Adiabatic growth in the isothermal sphere for masses = 0, 0.001, 0.01, 0.1, and 
1.0, with the mass increasing from bottom to top in the top two panels, and top to bottom 
in the last panel. Density is shown in the top frame, the radial velocity in the middle, and 
the anisotropy paramter /? in the bottom frame. 
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Fig. 2. — Adiabatic growth in the self-similar system with parameter 5 = 2 and masses 
iTLbh = 0, 0.001, 0.01, 0.1, and 1.0, with the mass increasing from bottom to top in the top 
two panels, and top to bottom in the last panel. Density is shown in the top frame, the 
radial velocity in the middle, and the anisotropy paramter (3 in the bottom frame. 
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Fig. 3. — Adiabatic growth in the self-similar system with parameter 5 = 3/4 and masses 
rribh = 0, 0.001, 0.01, 0.1, and 1.0, with the mass increasing from bottom to top in the top 
two panels, and top to bottom in the last panel. Density is shown in the top frame, the 
radial velocity in the middle, and the anisotropy paramter (3 in the bottom frame. 
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Fig. 4. — Adiabatic growth in the NFW system with masses m^/j = 0, 0.001, 0.01, 0.1, and 
1.0, with the mass increasing from bottom to top in the top two panels, and top to bottom 
in the last panel. Density is shown in the top frame, the radial velocity in the middle, and 
the anisotropy paramter (3 in the bottom frame. 
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Fig. 5. — — cr relationship in the three systems studied. Panel (a) shows the isothermal 
sphere, panel (b) is the self-similar system with S = 3/4, and panel (c) shows the NFW 
system. Shown is the relationship at three radii: R — 10~^ (solid line), R = 1 (dashed line), 
and R = 100 (dotted hne). 



